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ABSTRACT 

The  large  (V/m)  mesospheric  electric  fields  have  been  identified  as  a  possible  cause  of  VLF  phase 
perturbations.  These  fields  affect  the  fundamental  processes  that  govern  the  lower  D  region  parameters, 
primarily  the  electron  temperature  and  effective  collision  frequency.  The  main  ionospheric  parameter 
needed  to  calculate  VLF  phase  perturbations  is  the  low-frequency  electron  plasma  conductivity.  All  the 
electric  field  data  available  to  1990  were  collected  with  electric  field  sensors  on  board  more  than  50 
rockets  launched  over  approximately  30  years  in  the  USSR  and  the  U.S.A.,  which  were  insufficient  to 
address  VLF  phase  perturbations.  This  paper  discusses  the  progress  made  in  addressing  large  (V/m) 
mesospheric  electric  fields  between  60-  and  70-km  altitudes  since  1990.  It  focuses  on  achieving  the 
breakthrough,  the  development  of  a  radio  wave  technique  for  sensing  large  electric  fields  remotely  by 
using  MF  radar,  and  on  the  fact  that  the  electric  field  variability  leads  to  the  variability  of  ionospheric 
conduction  contours  by  a  few  kilometers  in  altitude.  The  statistical  analysis  of  the  large  mesospheric 
electric  field  data  acquired  in  the  60-  and  67 -km  altitude  region  in  Canada  and  Ukraine  suggests  that 
large  mesospheric  electric  fields  may  occur  during  about  70%  of  all  the  time.  However,  reasonable 
assessments  of  VLF  phase  perturbations  need  information  on  the  temporal  and  especially  spatial 
variability  of  conduction  contours,  which  remains  a  major  challenge  within  this  problem.  First,  the 
technique  developed  to  specify  electric  fields  requires  signal-to-noise  ratios  in  excess  of  a  factor  of  five, 
which  is  achieved  irregularly  with  the  MF  radars  used  at  present.  Second,  the  existing  MF  radars  do  not 
permit  the  observations  of  the  spatial  evolution  of  these  fields  at  all.  The  latter  problem  can  be  overcome 
by  developing  dedicated  radar.  Meanwhile,  co-located  VLF  phase  perturbation  measurements  and 
electric  field  observations  by  existing  MF  radars  may  be  combined  to  produce  a  pre-intermediate 
capability.  Eventually,  a  better  understanding  of  the  dynamics  and  mesospheric  and  ionospheric  D-region 
chemistry,  which  establish  conductivity  patterns,  will  require  the  combined  efforts  of  the  entire  scientific 
community. 
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ORGANIZATION 


1.0  INTRODUCTION 

The  mesosphere  is  an  active  electrodynamic  region  where  large  volt/meter  mesospheric  electric  fields  are 
generated  locally.  Much  of  the  evidence  for  and  against  large  mesospheric  electric  fields  exist  to  1990  is 
collected  in  the  reviews  by  Goldberg  [1-3].  All  the  data  to  1990  were  obtained  with  electric  field  sensors 
on  board  more  than  50  rockets  launched  over  about  30  years  in  the  USSR  and  the  U.S.A.,  and  they  were 
reported  in  tens  of  research  papers  [e.g.,  1-  4,  33  and  references  cited  therein]. 

Figure  1,  taken  from  the  paper  by  Maynard  et  al.  [4]),  shows  representative  vertical  electric  field  profiles 
from  rocket-borne  probes.  The  large  mesospheric  electric  fields  usually  occur  in  a  layer  of  order  10  km 
thick  above  about  60  km. 


Electric  Field  {Vim) 


Figure  1.  Electric  field  profiles  from  rocket-borne  symmetric  double  probes  at  Wallops  Island, 
Virginia,  on  31  July  1980  (after  Maynard  et  al.  [4]).  The  two  profiles  represent 
x  and  y  axis  sensors,  which  were  prepared  with  different  coatings. 


The  possible  physical  mechanisms  responsible  for  the  generation  of  large  mesospheric  electric  fields  have 
been  developed  in  a  few  papers  [5-8].  These  theoretical  models  provide  the  main  conditions  for  the 
occurrence  of  apparent  V/m  electric  fields  in  the  lower  mesosphere,  winds  and  heavy  ions. 

Seen  from  our  state  of  knowledge  in  1990,  a  remote  sensing  instrument  employing  a  radio-wave  technique 
was  critical  to  achieving  the  scientific  breakthrough  necessary  to  investigate  the  electrodynamics  of  the 
mesosphere. 

Since  1990,  we  have  developed  a  radio  wave  technique  for  remotely  sensing  large  volt/meter  electric 
fields,  which  are  intrinsic  to  the  mesosphere  [9  -  14]. 

Over  the  years,  the  MF  radars  in  Canada  and  Ukraine  have  accumulated  a  dataset  of  about  350  separate  5 
-10  min  intervals  of  measurements,  which  is  greater  then  the  entire  rocket  dataset  of  these  fields.  The  MF 
radar  dataset  has  yielded  knowledge  on  the  local  height  and  temporal  evolution  of  these  fields.  However, 
this  dataset  does  not  contain  at  all  information  about  the  three-dimensional  distribution  of  mesospheric 
electric  fields,  and  therefore  it  is  insufficient  to  model  VLF  phase  perturbations  thoroughly.  The 
determination  of  the  electric  field  spatial  distribution  poses  a  major  challenge  for  the  future. 
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The  next  section  will  present  the  MF  radar  technique  for  determining  large  mesospheric  electric  fields  and 
the  summary  of  its  capabilities.  Section  3  will  describe  the  disturbances  in  conductivity  profiles,  which 
emerge  from  modeling  studies.  Section  4  will  state  the  challenges  in  determining  VLF  phase  perturbations 
produced  by  large  mesospheric  electric  fields.  Section  5  will  indicate  the  immediate  benefits  of  co-located 
VLF  phase  perturbation  measurements  and  electric  field  observations  by  existing  MF  radars.  Finally, 
Section  6  will  summarize  the  achievement  and  problems  encountered  in  assessing  VLF  phase 
perturbations  produced  by  the  large  mesospheric  electric  fields. 


2.0  TECHNIQUE  FOR  DETERMINING  LARGE  MESOSPHERIC  ELECTRIC 
FIELDS 

2.1  Instrumentation 

The  MF  radar  at  the  Institute  of  Space  and  Atmospheric  Studies  (ISAS),  University  of  Saskatchewan, 
Canada  provided  polarimeter  data  for  the  61  -  67-km  altitude  range  using  a  20-ps  pulse  length  at  2.2  MHz 
during  1979  -  1982.  The  radar  is  a  single  frequency  system,  but  with  a  large  choice  of  antenna  arrays,  and 
is  well  described  in  recent  papers  [e.g.,  15].  It  is  used  mainly  for  wind  measurements,  but  has  also 
provided  studies  of  electron  densities  using  the  Differential  Absorption  Experiment  [16]. 

In  Ukraine,  the  measurements  were  made  with  the  Kharkiv  V.  Karazin  National  University  MF  radar  [17] 
the  specifications  for  which  are  as  follows:  operational  frequency  band  of  1.5-15  MHz,  16-element 
linearly  polarized  antenna  array  of  300  x  300  m2  physical  aperture  at /=  1.5^1. 5  MHz  and  of  60  x  60  m2  at 
/  =  4.5  -  5  MHz,  circularly  polarized  receiving  array  of  two-crossed  double  rhombus  antennas, 
polarization  switch  of  22  dB,  transmitter  peak  power  of  100  kW,  average  power  of  100  kW,  pulse  length 
of  20  ps  up  to  continuous  mode,  pulse  repetition  rate  of  1  -  100  s_1,  receiver  dynamic  range  of  86  dB,  IF 
bandwidth  of  60  kHz.  The  data  used  in  this  study  were  acquired  during  1978  through  1997  at  frequencies 
of /=  1.8  -  3.0  MHz  using  a  25-ps  pulse  length.  The  measurements  selected  for  this  study  were  made  in 
the  ionospheric  D  region  during  conditions  disturbed  only  with  respect  to  large  mesospheric  electric  fields 
and  quiet  with  respect  to  all  other  possible  disturbances,  when  the  signals  exceeded  the  noise  by  a  factor  of 
more  than  five.  The  observations  of  the  effective  electron  collision  frequency,  v,  were  made  by  applying 
the  “Differential  Absorption”  technique  of  [18]  at  the  altitudes  of  60-66  km. 

2.2  Basic  Relations 

The  basic  functional  relations  between  the  large  mesospheric  electric  field  features,  ionospheric 
characteristics,  and  scattered  signal  parameters  for  the  quasi-steady  case  are  given  by  [e.g.,  14,  19] 
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extraordinary  component  amplitudes  squared  and  averaged  over  the  sample; 

ordinary  component  amplitudes  squared  and  averaged  over  the  sample; 

total  absorption  of  the  ordinary  (+  subscript)and  extraordinary  (-  subscript) 
scattered  signal  components,  K_(z)  >  K+(z ); 

high  frequency  electron  conductivity  for  the  ordinary  and  extraordinary 
components,  respectively; 

2  if, 

sounding  frequency; 

2  4l, 

electron  gyrofrequency  longitudinal  component; 

1.35  MHz  for  middle  latitudes; 

coefficient  that  represents  the  kinetic  effects  in  the  high  frequency  permittivity 
[e.g.,  20]; 

coefficient  that  represents  the  kinetic  effects  in  the  high  frequency  conductivity 
[e.g.,  20]; 

electron  charge; 
electron  mass; 
speed  of  light; 

lower  ionospheric  boundary  altitude  that  is  assumed  to  be  equal  to  z0  =  60  km 
for  the  daytime  and  z0  =  75  -  80  km  for  the  nighttime,  in  most  of  the  cases; 
ion  production  rate; 

effective  rate  at  which  negative  ions  are  destroyed  by  electron  detachment; 
electron  number  density; 

JST/N 

negative  ion  number  density; 

effective  rate  at  which  the  negative  ions  are  formed  by  attachment  of  electrons 

to  neutral  constituents; 

effective  rate  of  electron-ion  recombination; 

effective  rate  of  ion-ion  recombination; 

mean  energy  imparted  to  an  electron  by  an  external  heating  source,  e.g., 
external  electric  field; 

Boltzmann's  constant; 
electron  temperature; 
neutral  species  temperature; 

fractional  loss  of  energy  per  electron  collision  with  a  heavy  particle; 
effective  electron-neutral  collision  frequency; 

density  of  the  current  driven  by  an  external  atmospheric  current  source; 
low-frequency  conductivity  of  the  ionospheric  Z)-region  plasma; 
quasi-steady  vertical  mesospheric  electric  field  intensity. 
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Here,  Equation  (1)  represents  the  relation  between  the  disturbed  Z)-region  parameters  and  the  ratio  of  the 
squared  relative  amplitudes  of  the  ordinary  and  extraordinary  components  of  the  scattered  signals  in  the 
quasi-longitudinal  approximation.  Equations  (4)  and  (5)  are  the  continuity  equations  for  the  electrons  and 
ions,  respectively,  (6)  is  the  energy  equation,  and  (7)  is  nonlinear  Ohm's  law  for  the  quasi-steady  large 
mesospheric  electric  fields.  In  writing  Equations  (4)  -  (6),  it  has  been  assumed  that  the  weakly  ionized 
ionospheric  plasma  is  quasi-neutral,  the  positive  and  negative  ion  temperatures  are  equal  to  the  neutral 
constituent  temperature,  and  the  effects  of  transport  processes  on  local  disturbances  can  be  neglected  [e.g., 
14].  In  the  D  region 

Q,=j,.E  =  jll<je.  (8) 


Also,  the  following  are  taken  into  account  [e.g.,  20,  21]: 
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where  KG  =  1.42  [20],  Nn  is  the  number  density  of  neutral  particles,  N(02)  is  the  number  density  of 
molecular  oxygen  in  cm  3, 7V(N2)  is  the  number  density  of  molecular  nitrogen  in  cm-3,  Te  and  Tn  are  in  K, 
va  in  s  ’,  ar  in  cm3  s  ',  the  subscript  0  is  used  to  denote  the  magnitude  of  the  plasma  parameters  in  the 
absence  of  large  mesospheric  electric  fields. 


When  the  differential  absorption  of  the  two  magnetoionic  components  is  neglected  in  the  first,  the  lowest 
sample  of  signals  scattered  from  the  layer  at  the  altitude  zi,  which  is  usually  true  for  z\  <  66  -  69  km  for 
the  daytime  conditions  [e.g.,  18],  the  following  relation  is  derived  from  (1)  for  determining  the  disturbed 
ve(z\)  value  from  the  measured  value  R(z\)  [14,  22,  13]: 
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Then,  from  (6)  -  (1 1)  it  is  easy  to  obtain  a  relation  for  the  electric  field  intensity  E(z\)  [14,  22] 
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where  Te0(z\)  and  veo(z\)  are  related  by  (10). 
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ORGANIZATION 


The  determination  of  the  current  density  j  =  constant  by  using  (4)  -  (16)  is  a  crucial  step  before  the 
transition  to  the  second  layer  from  which  the  signals  have  been  scattered.  This  condition  is  apparently 
satisfied  in  the  60  to  70  -  75-km  altitude  range  (see  [12]  for  a  discussion  of  this  subject).  The  disturbed 
value  of  N(z\)  is  given  by 


Nizi )  =  q'/2  (zi ){(! +HQ  ))(«,■  (#i )+ a(6>  )«,.  (z, ))} 


(17) 


which  allows  the  specification  of  je  by  (7),  (9).  Here,  9\  =  Te{z\ )  /  Te0(zi)  =  ( ve(z\ )  /  ve0(zi))615.  Then,  from 
(6)  -  (9)  the  N(z2 )  in  the  second  layer  is  given  by 


N(z2)  =  je\2.lxk-Te0(z2)S(O2)(92-l) 
m 


(18) 


Equating  (18)  and  the  expression  for  N(z2)  derived  from  (1)  gives  the  following  equation  for  ve(z2) 
allowing  for  the  differential  absorption  of  the  scattered  signals: 
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where  A z  is  the  height  increment.  Since  both  the  right-hand  and  left-hand  sides  of  (19)  are  equal  to  N(z2), 
then  the  N(z2)  can  be  easily  determined  provided  the  value  of  v(z2)  has  been  established. 

We  obtain  the  following  equation  for  v(z3)  in  the  third  layer: 

je\2.U-Te0{z2)5(92)(d2-\) 

1  m 

r,  +2  2  (  Kl  UJ~UJ>  +v2 (zOk2  ^  x 

\{co  +  coL)  +ve(z2)\y  O  £{ve(z3))  y  , 


ln< 


f  \ 

\ 2  __9  CO  -CO 


CO  -  CO, 


[(co-coL)2+v2e(z3)]  ^  +  (0l) 


2k2 


c  \ 

CO  +  CO 


XzX 


+v2Az3)Kl 


c  \ 

cox- CO 


XzA 


>-4  (k_(z2)-K+(z2)) 


%ne 


;(z3)Az 


CO-COr 


•(Z3) 


y 

V  ^  \~ijj 


co  +  co. 


.(z3). 


(w-0)L)  +  v2  (z3)  (co  +  c vL)  +V2(zl) 


(20) 
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where 


4  {K_(z2)-K+(z2j)  =  ln< 


R~'(z2) 


[(®  +  ®z)2+v2(z2)] 

S 

i 

f  \ 

CO  -  CO, 

U±)J 

r  \\ 

co  —  co L 

U(z2)J 

[(a)-aL)2  +v2e(z2jj 

( (o  +  a>L)2  K) 

f  \ 

CO  +  CD, 

U±)J 

+v2(*2K 

f  \ 

CO  +  CO, 

U±)J. 

is  obtained  from  (1)  -  (3),  (19).  The  N{z3)  is  obtained  from  (20)  as  above  by  using  the  value  of  ve{z3) 
determined  from  (20). 

Further,  the  procedure  is  consecutively  repeated  until  the  relative  perturbations  in  the  effective  electron 
temperature  and  collision  frequency  become  much  less  than  unity.  Then,  the  N  in  the  highest  layers  that 
are  practically  undisturbed  is  determined  from  (1),  (2)  by  applying  the  classical  differential  absorption 
technique  and  specifying  the  value  of  ve0  [18]. 

It  is  important  to  note  that  this  technique  permits  not  only  simultaneous  measurements  of  the  height 
profiles  of  large  mesospheric  electric  field  intensities  and  the  profiles  of  associated  disturbances  in  the 
effective  collision  frequency,  the  electron  temperature  and  density,  but  also  simultaneous  estimates  of 
disturbances  in  the  D-region  basic  parameters  va/vao,  8/80,  ar/ar0,  AA0,  a/aeo,  and  in  N~  /  and  N+  /  A0+ . 
Thus,  the  developed  technique  allows  remote  sensing  of  a  cluster  of  the  D-region  parameters  disturbed  by 
the  large  mesospheric  electric  fields.  Unfortunately,  an  essential  deficiency  of  this  technique  is  the  fact 
that  it  is  applicable  only  when  the  strength  of  the  signals  scattered  from  the  lower  part  of  the  D  region,  60 
-  70  km,  is  high  enough,  which  does  not  always  occur  when  the  existing  MF  radars  are  used. 

2.3  Technique 

This  section  illustrates  how  the  technique  works  in  detail.  The  data  selected  for  this  analysis  were 
collected  on  November  24,  1984  with  the  MF  radar  at  the  Kharkiv  V.  Karazin  National  University 

o  o 

Radiophysical  Observatory  (49  38?N,  36  20?E)  [e.g.,  17].  The  radar  specifications  in  this  experiment  were 
as  follows:  2.3-MHz  sounding  frequency,  25-ps  pulse  length,  and  pulse  repetition  rate  of  1  per  second.  We 
have  chosen  for  the  analysis  the  5-min  interval  09:54  -  09:59  LT,  which  exhibits  no  time  trend  in  the 
intensity  of  the  signals  scattered  from  the  66  -  84-km  altitude  and  the  signal-to-noise  ratio  exceeding  8. 

After  testing  the  input  statistical  series  of  the  squared  noise  amplitudes  and  of  the  squared  signal  plus 
noise  amplitudes  for  homogeneity  and  further  subtracting  the  noise  from  the  signal  plus  noise,  the  height 
dependence  of  R(z)  defined  in  (1)  is  formed,  and  the  corresponding  sample  variances  S2(R)  are  estimated 
[e.g.,  23]: 


s2(R) 


WV(£)+£(£)  4g)4g)| 

At  1  *  Al  ±  Al-  Al  J 


(21) 


where  S2^A2_^,  represent  the  sample  variance  of  the  extraordinary  and  ordinary  scattered  signal 

components,  respectively,  p±  is  the  sample  correlation  coefficient  between  the  squared  extraordinary  and 
ordinary  scattered  signal  components.  Equation  (21)  yields  the  99  percent  confidence  intervals  for  7?(z), 
which  are  lying  within  the  range  ±2.4%  for  7?(66  km)  and  ±16%  for  R( 78  km). 


Then,  the  dependence  R{z)  obtained  experimentally  is  used  for  simultaneously  determining  the  height 
dependences  of  the  basic  parameters  of  the  disturbed  D  region  and  the  disturbing  large  mesospheric 
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electric  field  characteristics  shown  in  Table  1  where  the  0.99  confidence  intervals  are  indicated  in 
parentheses. 

First,  the  ve(66  km)  is  determined  by  making  use  of  (14)  and  veo(66  km)  =  1.68  x  107  s”1  that  corresponds 
to  an  undisturbed  electron  temperature  of  Te0  =  240  K  in  the  isothermal  ionospheric  plasma;  the  model 
parameter  values  of  the  undisturbed  medium  are  borrowed  from  [20,  21,  24  -  30].  Then,  given  the  ve,  ve0, 
the  values  of  E ,  TJTe 0,  S/So,  a/cco,  v/vo,  AJAo  for  z  =  66  km  are  obtained  from  (10)  -  (13),  (15).  Taking  into 
account  the  high  stability  of  the  main  ionizing  radiations,  the  cosmic  rays,  at  66  km  under  quiet 
geomagnetic  conditions  [e.g.,  24  -  26,  30,  31]),  the  electron  number  density  7V0(60  km)  can  be  set  equal  to 
60  cm-3  for  qt  =  0.05  cm-3  s”1  (see  also  [12]).  Then,  from  (4),  (5),  (12),  (13)  when  aro  =  6.7  x  1(T6  cm-3  s”1, 
at  6.8  x  10”8  cm-3  s”1,  va0  =  0.7  s”1,  vd  =  0.7  s”1,  Ao  va0/ vd  =  1.0,  we  obtain  N(66  km),  and  from  (7),  (9) 
we  also  find  the  height-independent  current  density  je  =  3.58  x  1(T8  A/m2  induced  by  the  mesospheric 
source  at  z  >  66  km. 

The  disturbed  values  of  the  ionospheric  plasma  parameters  at  higher  altitudes  are  determined  in  the  same 
way  as  above.  The  initial  values  are  assumed  to  be  equal  to 

ve0  (69  km)  =  1.1  x  107  s"1,  Te0  (69  km)  =  230  K, 

ve0  (72  km)  =  6.8  x  106  s"1,  Te0  (72  km)  =  210  K, 

ve0  (75  km)  =  4.2  x  106  s"1,  Te0  (75  km)  =  200  K. 

Table  1  shows  that  the  relative  measurement  errors  in  the  ionospheric  parameters  derived  from  the  given 
experiment  data  vary  within  the  limits  of  the  order  of  one  percent  at  z  =  66  km  to  23%  at  z  =  75  km,  as 
determined  for  the  99  percent  confidence  intervals.  However,  there  exist  sources  of  bias  errors  that  require 
analysis  in  each  specific  event,  e.g.,  in  this  experiment,  an  additional  error  of  2%  is  due  to  the  neglect  of 
the  differential  absorption  at  z  <  66  km.  In  general,  the  assumed  values  of  ve0  and  associated  with  them  Te0 
(see  (10))  also  exhibit  variations  due  to  seasonal,  diurnal,  and  regional  dependences  of  atmospheric 
parameters,  as  well  as  the  atmospheric  disturbances  from  acoustic-gravity  waves  [e.g.,  32].  The  errors 
caused  by  radar  imperfections  may  not  also  be  excluded  altogether  in  determining  ve  (e.g.,  inaccuracy  in 
altitude  determination).  The  control  and  taking  account  of  all  possible  errors  is  very  difficult.  However, 
Equation  (1)  shows  that  each  specific  sample  value  of  R  has  an  associated  interval  of  ve0  that  corresponds 
to  the  positive  values  of  A  in  (1).  Therefore,  the  physical  requirement  N  >  0  actually  imposes  a  restriction 
on  the  maximum  possible  total  bias  in  determining  or  assuming  the  value  of  veo.  For  example,  this  kind  of 
error  in  this  experiment  is  not  more  than  16%,  which  corresponds  to  N  =  0  at  z  =  69  km.  The  analysis  of 
the  Kharkiv  MF  radar  data  has  revealed  that  this  kind  of  error  attained  a  maximum  value  of  38%  over  the 
interval  1979  -  1994  only  once  (probability  of  0.005)  [13,  22].  The  mean  was  approximately  equal  to 
19%.  The  inverse  proposition  is  also  valid:  a  limited  interval  of  R  ±  AR  corresponds  to  the  existing  interval 
of  veo  ±  Ave0,  which  limits  the  maximum  possible  bias  in  determining  R ;  in  the  given  experiment  the 
maximum  possible  value  of  AR  /  R  is  equal  to  approximately  16%.  Thus,  in  a  specific  experiment,  there 
always  exist  physical  constraints  for  restricting  and  estimating  the  maximum  possible  bias.  For  example, 
in  this  experiment,  the  bias  does  not  exceed  30%  at  z  <  69  km  and  50%  in  the  interval  69  km  <  z  <  75  km, 
which  is  in  accordance  with  the  results  of  Gokov  and  Martynenko  [22],  Martynenko  et  al ,  [13],  and 
Martynenko ,  Rozumenko ,  and  Tyrnov  [12]. 

As  a  whole,  Table  1  shows  that  large  mesospheric  electric  fields  cause  significant,  by  a  factor  of  up  to  two 
times  and  more,  disturbances  virtually  in  all  basic  parameters  of  the  lower  Z)-region  plasma,  except  for  the 
electron  number  density.  The  latter  is  explained  by  a  mutual  balance  between  the  processes  resulting  in  a 
decrease  in  N  due  to  the  attachment  of  electrons  to  air  molecules  and  the  processes  resulting  in  an  increase 
in  N  due  to  a  decrease  in  the  effective  rate  of  electron-ion  recombination  under  the  influence  of  large 
mesospheric  electric  fields  (see  (4)  -  (7),  (12),  (13)). 
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Table  1  Parameters  of  the  D  region  disturbed  by  large  mesospheric  electric  fields. 


z  (km) 


Ve  (S'*) 


E  (V/m) 


Ve  /  Ve0 


Te/T, 


&  1  e0 


5/So 


N  (cm  ) 


Va  /  Va0 


ar/a. 


T/Ur0 


N/N0 


N-/N~ 


N+  /  Nq 


ae/a 


e/'J  e0 


66 


3.15 


+0.07 

-0.07 


xlO7 


0.46 


+  0.02 
-0.02 


1.88 


+0.05 

V-0.05y 


2.13 


+0.06 

-0.06 


0.47 


+0.01 

-0.01 


60 


^+0 

V-2y 


1.95 


+0.06^ 

-0.06y 


0.69 


^+0.01 

V-0-01y 


1.95 


+0.06 

-0.06 


1.00 


^+0.01 

v-0.03y 


1.95 


+0.01 

-0.04 


1.50 


+0.01 

-0.04 


0.53 


+0.02 

-0.03 


69 


2.51 


+0.22 

v-0.29y 


xlO7 


0.40 


+0.04 

-0.07 


2.29 


+0.19 

-0.27 


2.70 


+0.28 

v-0.38y 


0.37 


+0.06 

-0.03 


60 


^+5  ^ 
V  12y 


2.15 


f +0.03 
^ — 0.05  j 


0.61 


+0.05 

V-0.03y 


2.15 


+0.03 

-0.05 


1.11 


+0.09 

-0.22 


2.39 


+0.13 

V-0.45y 


1.47 


+0.11 
v-0.28  j 


0.48 


+0.01 

-0.04 


72 


6.90 


+0.10 

-0.10 


xlO6 


0.02 


+0.01 

-0.01 


1.01 


+0.01 

-0.01 


1.01 


+0.01 

-0.01 


0.99 


+0.01 

-0.01 


350 


^+50^ 

V 60y 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


0.99 


+0.01 

v-O-Oly 


75 


4.2x10 


1.0 


1.0 


1.0 


400 


^+90^ 

V 80y 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


The  technique  developed  for  remotely  sensing  large  mesospheric  electric  fields  can  briefly  be  summarized 
as  follows: 

(i)  The  ve{z\)  is  determined  by  using  (14)  in  the  first  lowest  layer  at  z\  <  66  -  69  km  from  which  the 
signals  have  been  received, 

(ii)  The  Te(zi)/Te0(zi)  is  specified  by  using  the  ve{z\)  and  (10), 

(iii)  The  E(zx)  is  determined  by  using  (15)  or  (16)  depending  on  the  value  of  Te(z\ )  /  Te0(zi), 

(iv)  The  N(z{)  is  calculated  by  using  (17), 

(v)  The  current  density  je  =  constant  for  the  source  of  large  mesospheric  electric  fields  is  determined  from 
(7)  and  (9), 
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(vi)  The  ve(z2)  is  determined  for  the  second  layer  of  the  scattered  signals  from  (19)  and  then  the  N(z2)  is 
determined, 

(vii)  The  E(z2)  is  determined  from  (7)  and  (9), 

(viii)  The  ve(z3)  is  determined  in  the  third  layer  of  the  scattered  signals  from  (20)  and  then  the  N(z2)  is 
determined, 

(ix)  The  E(z2)  is  determined  from  (7)  and  (9), 

(x)  Further,  the  procedure  is  consecutively  repeated  by  using  equation  (20)  and  relations  (7),  (9)  until  the 
relative  disturbances  in  the  electron  temperature  and  effective  collision  frequency  become  much  less  than 
unity;  at  higher  altitudes  the  medium  can  be  considered  undisturbed  and  the  values  of  E  small,  i.e.,  close  to 
zero. 

Hence,  the  set  of  theoretical  relations  (1)  -  (17)  provides  the  framework  for  modeling  studies  of  how  large 
mesospheric  electric  fields  affect  the  ionospheric  D-region  parameters.  The  disturbances  in  the  electron 
temperature  and  effective  collision  frequency  (see  Equations  (10),  (15),  (16))  are  the  primary  cause  of 
disturbances  in  other  parameters.  In  particular,  Equations  (1),  (14)  describe  disturbances  in  the  ratio  of  the 
squared  amplitudes  of  the  ordinary  and  extraordinary  components  of  the  scattered  signals.  Equation  (2) 
governs  variations  in  the  total  absorption  of  MF  radio  signals.  Equation  (3)  permits  the  determination  of 
disturbances  in  the  HF  conductivity  of  the  ionospheric  plasma.  Equation  (7)  relates  the  large  mesospheric 
electric  fields  and  the  low-frequency  electron  conductivity  of  the  plasma.  Equation  (8)  takes  account  of  the 
large  mesospheric  electric  field  energy  losses  via  Joule  heating.  Equation  (10)  provides  the  relationship 
between  disturbances  in  the  electron  temperature  and  the  effective  collision  frequency.  Equation  (11) 
establishes  disturbances  in  the  fractional  loss  of  energy  per  electron  collision  with  a  heavy  particle. 
Equation  (12)  is  used  to  calculate  the  effective  rate  at  which  the  negative  ions  are  formed  by  the 
attachment  of  electrons  to  neutral  constituents.  Equation  (13)  shows  disturbances  in  the  effective  rate  of 
electron-ion  recombination.  Equation  (17)  defines  explicitly  disturbances  in  the  electron  number  density. 

Generally,  the  theoretical  model  outlined  above  forms  the  basis  for  the  technique  of  Martynenko  [19]  for 
clustered-technique  remote  sensing  of  processes  coupling  the  electrically  active  mesosphere  with  the 
ionospheric  Z)-region  plasma,  the  accuracy  of  which  is  not  inferior  to  the  corresponding  in  situ  rocket 
techniques  [e.g.,  33]  but  offers  considerable  cost  benefits.  A  significant  extension  to  the  earlier  technique 
of  Gokov  and  Martynenko  [22]  and  Martynenko  et  al.  [13]  has  been  achieved  by  using  the  model 
representation  of  the  source  of  large  mesospheric  electric  fields  as  the  current  source  [12].  It  includes  the 
extension  of  the  altitude  range  and  permits  the  simultaneous  measurements  of  large  mesospheric  electric 
fields  and  charged  particle  number  densities.  The  minimum  total  error  of  the  suggested  method  is 
approximately  20%  and  the  maximum  total  error  can  attain  a  value  of  40  -  50%,  while  the  particular  error 
magnitudes  need  to  be  specified  in  each  particular  experiment.  The  above-mentioned  errors  can  be 
reduced  by  roughly  2  times  if  independent  measurements  of  the  neutral  temperatures  in  the  mesosphere 
were  provided,  e.g.,  by  the  lidar  techniques  [34]. 

Thus,  there  are  good  grounds  for  the  conclusion  that  further  development  and  implementation  of  the 
remote  sensing  techniques  will  reduce  an  acute  shortage  of  data  on  the  electrodynamic  processes  acting  in 
the  electrically  active  mesosphere  and  on  the  mesosphere's  coupling  to  other  atmospheric  and  ionospheric 
regions.  These  techniques  will  permit  simultaneous  measurements  of  large  mesospheric  electric  fields  and 
disturbed  lower  ionospheric  parameters. 

2.4  Large  Mesospheric  Electric  Field  Measurements 

This  section  briefly  summarizes  the  main  results  inferred  thus  far  from  MF  radar  electric  field 
measurements  taken  in  Canada  and  Ukraine  [9,  10,  11].  Below,  we  characterize  the  database,  illustrate 
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temporal  variability,  and  present  the  statistical  analysis  of  the  large  mesospheric  electric  fields,  electron 
temperatures,  and  the  effective  collision  frequencies. 


2.4.1  Database 

The  estimates  of  VLF  phase  variations  require  the  characterization  of  the  temporal  and  3 -dimensional 
variability  of  the  ionospheric  conduction  contours.  However,  the  ionospheric  observations  have  so  far 
been  made  with  monostatic  MF  radars  with  antennas  forming  wide  fixed  vertical  beams  and  thus  capable 
of  measuring  only  a  one-dimensional  distribution  of  ionospheric  parameters.  The  horizontal  variability 
remains  generally  unknown.  Consequently,  the  available  MF  radar  database  is  insufficient  to  characterize 
the  lower  ionospheric  boundary  in  detail  in  order  to  compute  VLF  phase  variations. 

The  database  of  MF  radar  measurements  is  developed  in  a  series  of  campaigns  carried  out  in  Ukraine 
between  1978  -  1994  and  in  Canada  between  1979  -  1982.  The  large  mesospheric  electric  field  and  the 
effective  electron  collision  frequency,  v;  data  were  taken  from  the  60  -  66  km  altitude  range  at  Kharkiv  V. 
Karazin  National  University  and  from  the  61  -  67  km  range  at  the  Institute  of  Space  and  Atmospheric 
Studies  (ISAS),  University  of  Saskatchewan,  Canada.  The  Ukrainian  MF  radar  operated  at  1.8  -  3.0-MHz, 
and  acquired  n  =  185  measurements  of  5  -  10  min  in  duration  each  [13,  22].  The  Canadian  MF  radar 
operated  at  2.2  MHz,  and  collected  n=  \10  measurements  of  10  min  in  duration  each.  Thus,  the  total  MF 
radar  database  contains  385  MF  radar  measurements,  which  exceeds  by  a  factor  of  a  few  times  the  total 
rocket  database. 

2.4.2  Temporal  Variability 

The  variability  of  large  mesospheric  electric  fields  displays  a  wide  range  of  time  scales,  from  an  order  of 
one  minute  (as  shown  in  Figure  1  in  [10])  to  an  order  of  an  hour,  as  shown  in  Figure  2  for  Ukraine 
(probing  frequency  of '/=  2.3  MHz)  and  in  Figure  3  for  Canada  (f=  2.2  MHz). 

However,  the  database  consists  of  separate  5-10  min  intervals  of  measurements  collected  over  the  years 
and  the  temporal  variability  they  provide  is  not  continuous.  The  continuity  of  measurements  may  be 
improved  by  either  upgrading  the  existing  MF  radars  or  constructing  new  radars  in  order  to  increase  the 
signal-to-noise  ratio. 

2.4.3  The  Distribution  of  the  Effective  Electron  Collision  Frequencies 

The  effective  electron  collision  frequency,  v,  is  the  primary  plasma  parameter  that  is  used  for  deriving 
both  the  electric  fields  and  the  low-frequency  conductivities.  Preliminary  data  analysis  includes  the 
estimate  of  the  histogram  showing  the  distribution  of  v  /vm  where  vm  is  a  model  value  of  v  in  the  absence 
of  large  mesospheric  electric  fields  at  heights  from  which  the  scattered  signals  were  received.  Allowing 
for  the  kinetic  effects,  as  Gokov  and  Martynenko  [22]  and  Martynenko  et  al.  [14]  assumed,  vm  =  3.75  x  107 
s”1  at  z  =  60  km,  vm  =  3.32  x  107  s”1  at  z  =  61  km,  vm  =  2.55  x  107  s”1  at  z  =  63  km,  vm  =  2.21  x  107  s”1  at  z 
=  64  km,  vm  =  1.68  x  107  s”1  at  z  =  66  km,  and  vm  =  1.47  x  107  s”1  at  z  =  67  km.  The  histograms  showing 
the  distribution  of  v/vm  values  are  presented  in  [10]  (Figure  3  from  the  Ukrainian  data  and  Figure  4  from 
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Figure  2.  A  representative  time  dependence  of  large 
mesospheric  electric  fields  obtained  in  Ukraine. 
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Figure  3.  A  representative  time  dependence  of  large 
mesospheric  electric  fields  obtained  in  Canada. 


the  Canadian  data).  It  is  obvious  that  the  v/vm  sample  in  the  absence  of  large  mesospheric  electric  fields 
should  have  a  Gaussian  distribution  with  a  mean  of  M  [  W  vm]  =  1.  This  provides  the  reason  for  suggesting 
that  the  values  W  vm  <  1  are  associated  with  cases  when  large  mesospheric  electric  fields  are  absent,  and 
these  data  are  excluded  from  the  database  of  the  W  vm  measurements  that  are  used  for  determining  large 
mesospheric  electric  field  effects.  Furthermore,  the  symmetry  of  the  Gaussian  distribution  of  W  vm  about 
v  /  vm  =  1  in  the  absence  of  large  mesospheric  electric  fields  provides  another  criterion  for  excluding  the 
undisturbed  component  with  v!  vm>\  from  the  analysis. 
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2.4.4  The  Distribution,  wE ,  of  the  Large  Mesospheric  Electric  Fields 

The  histograms  for  the  distribution,  wE,  of  the  large  mesospheric  electric  fields,  E  are  shown  in  [10] 
(Figure  5  for  the  Ukrainian  and  Figure  6  for  Canada).  There,  wE  =  nt  Kyi  AE ),  i  is  a  cell  number,  i  = 
1,. .  .,13,  Yii  is  the  number  of  samples  of  E  that  lie  within  the  (i-  1)  AE  <E  <{i  AE)  cell,  A  E  =  0.5  V/m  is 
the  width  of  each  cell,  n  =  \29  (Ukraine)  and  n  =  120  (Canada)  are  the  sample  sizes.  Within  the  0.99 
confidence  interval,  the  histograms  exhibit  the  possibility  of  dividing  them  into  two  parts:  that  constituting 
the  main  body  of  0  <  E  <  2.5  V/m,  n=  \29  for  the  Ukrainian  data  and  n  =  99  for  the  Canadian  data,  and 
that  constituting  the  tail  E  >  2.5  V/m,  n  =  10  for  the  Ukrainian  data  and  n  =  21  for  the  Canadian  data.  Also, 
within  the  0.99  confidence  interval  estimated  by  making  use  of  Pearson's  test,  the  main  body  corresponds 
to  a  one-parameter  Rayleigh  probability  density  functions  as  given  by 


HE) 


(22) 


where  a  =  Em  is  the  most  probable  value  of  E ,  MX[E ]  =  ( n /  2) 1/2  a  is  the  mean,  M2[E\  =  2  o2  is  the  second 
ordinary  moment  for  a  Rayleigh  set,  and  D  [E]  =  (2  -  n  /  2)  o2  is  the  variance.  Within  the  same  0.99 
confidence  interval,  the  Rayleigh  part  of  the  histogram  in  [10]  (Figure  5)  provides  an  estimate  of  M\  [E]  = 
0.72  ±0.11  V/m  where  the  sample  mean  <  ER>  =  0.72  V/m,  which  corresponds  to  cr=  0.57  V/m.  The 
corresponding  histogram  wE  for  large  mesospheric  electric  field  intensities  constructed  by  using  n  =  \29 
samples  and  the  theoretical  Rayleigh  probability  density  function  are  presented  for  Ukraine  in  [10]  (Figure 
7).  The  similar  histogram  for  Canada  is  presented  in  [10]  (Figure  8)  where,  MX[E]  =  0.89  ±  0.12  V/m,  <  ER 
>  =  0.89  V/m,  and  a=  0.71  V/m. 

Hence,  at  least  two  mechanisms  for  generating  large  mesospheric  electric  fields  should  exist.  The 
Rayleigh  distribution  of  E  can  be  formed  as  a  result  of  the  summation  of  the  random  fields  from  a  large 
number  of  primary  small-scale  mesospheric  generators.  The  possible  mesospheric  processes  resulting  in 
such  small-scale  active  elements  have  been  discussed,  e.g.,  by  Goldberg  [3]  and  Polyakov  et  al.  [6].  The 
processed  data  have  shown  that  the  probability  of  occurrence  of  such  an  integral  Rayleigh  mesospheric 
generator  is  equal  to  approximately  70%  for  Ukraine  and  approximately  58%  for  Canada,  and  the 
probability  of  the  lack  of  large  mesospheric  electric  fields  is  about  25%  for  Ukraine  and  approximately 
30%  for  Canada. 


Unfortunately,  the  number  of  observations  occurring  in  the  sample  interval  E  >  2.5  V/m  is  equal  to  n  =  10 
with  the  probability  of  occurrence  of  approximately  5%  for  Ukraine  (see  Figure  5  in  [10])  and  n  =  21  with 
the  probability  of  occurrence  of  approximately  12%  for  Canada  (see  Figure  6  in  [10]).  This  dataset  is  too 
small  to  draw  statistical  inferences  from,  and  these  data  are  barely  adequate  for  characterizing  a  mean  of  < 
E  >  =  (4.3  ±  1.3)  for  Ukraine  and  <E>  =  (4.4  ±  0.4)  for  Canada  within  the  0.90  confidence  interval. 

Generally,  the  upper  limit  for  E  is  determined  by  an  atmospheric  breakdown  threshold  of  £)  =  218  x  (p  / 
Po)  kV/m  (here  p  is  the  atmospheric  pressure  at  the  altitude  z,  p0  is  the  atmospheric  pressure  at  sea  level 
[e.g.,  35]).  For  example,  Et  =  50  -  10  V/m  within  the  altitude  range  of  z  =  60  -  70  km,  respectively.  The 
fields  exceeding  this  threshold  are  associated  with  red  sprite,  blue  jet,  and  elf  phenomena  in  the  middle 
atmosphere. 


2.4.5  Seasonal  Dependencies  in  the  Statistical  Parameters  of  Large  Mesospheric  Electric  Fields 

To  detect  a  possible  seasonal  dependencies  in  the  statistical  parameters  of  large  mesospheric  electric 
fields,  all  the  data  in  the  database  were  arbitrary  divided  into  two  subsets:  "winter"  (September  24  -  March 
23)  and  "summer"  (March  24  -  September  23).  Then  the  Rayleigh  components  in  the  E  distribution  were 
constructed,  with  n  =  69  for  the  winter  and  n  =  60  for  the  summer  for  Ukraine  and  n  =  49  for  the  winter 
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and  n  —  50  for  the  summer  for  Canada.  The  corresponding  histograms  showing  the  distributions  of  wE  and 
the  theoretical  Rayleigh  distributions  are  presented  in  [10]  (Figure  11  and  Figure  12  for  Ukraine  and 
Figure  13  and  Figure  14  for  Canada).  Within  the  0.99  confidence  interval,  in  accordance  with  the 
Pearson's  test,  the  histograms  of  wE  correspond  to  Rayleigh  probability  density  functions  with  <  ER>  = 
0.70  V/m  and  <j=  0.56  V/m  for  the  winter,  and  <ER>  =  0.75  V/m  and  <j=  0.60  V/m  for  the  summer  for 
Ukraine.  For  Canada,  <ER>  =  0.91  V/m  and  cr=  0.73  V/m  for  the  winter,  and  <ER>  =  0.86  V/m  and  cr= 
0.69  V/m  for  the  summer.  It  can  be  seen  that  that  the  seasonal  differences  have  turned  out  to  be 
statistically  insignificant.  The  E  >  2.5  V/m  values  were  also  occasionally  observed  during  both  the  winter 
(?z  =  6  for  Ukraine  and  n  =  1 1  for  Canada)  and  the  summer  (ft  =  4  for  Ukraine  and  n  =  10  for  Canada). 

Generally,  this  may  indicate  that  the  mean  local  seasonal  variations  in  mesospheric  parameters,  for 
example,  due  to  the  mean  local  thunderstorm  activity  do  not  exert  a  noticeable  effect  at  least  on  the  mean 
performance  of  the  Rayleigh  generator  of  large  mesospheric  electric  fields. 

2.4.6  Diurnal  Dependence  in  the  Large  Mesospheric  Electric  Field  Statistics 

An  attempt  has  been  made  to  reveal  a  diurnal  dependence  in  the  large  mesospheric  electric  field  statistics 
in  [10].  The  histograms  for  the  Canadian  data  (Figures  15  and  16  in  [10])  that  show  the  distribution  of  wE 
are  fitted  with  theoretical  Rayleigh  distributions  for  the  day  (ft  =  72,  <  ER  >  =  0.91  V/m,  cr=  0.73  V/m) 
and  for  the  night  (ft  =  27,  <  ER  >  =  0.86  V/m,  cr=  0.69  V/m),  respectively.  A  comparison  of  these  results 
indicates  the  absence  of  a  noticeable  diurnal  dependence  of  the  distribution  function  of  large  mesospheric 
electric  field  values. 

2.4.7  Distribution  Functions  for  the  Relative  Disturbances  in  the  Effective  Electron  Collision 
Frequencies  and  in  the  Electron  Temperatures 

The  distributions  of  the  effective  electron  collision  frequencies  and  the  electron  temperatures  are  reflected 
in  the  distributions  of  disturbances  in  the  low-frequency  electron  conductivity  of  the  plasma. 

For  the  Rayleigh  distribution  of  large  mesospheric  electric  fields,  Equation  (22),  we  have  derived  the 
theoretical  distribution  functions  f{rj)  and  f{9)  for  the  relative  disturbances  in  the  effective  electron 
collision  frequency  r/  =  ve/  ve0  and  in  the  electron  temperature  9  =  Tet  Te0  by  making  use  of  the 
deterministic  functional  dependences  in  Equations  (10),  (16),  as  given  by 


f{0)  =  ^r(5^2/3  -  2d~m )exp  j-^2/3  (0  - !)}  >  (24) 

where  Sl  (z)  =  (kmS0Te0v^/(0. 97 e2 )  ,J{rj)  =  0  for  rj  =  1  ,f(9)  =  0  for  9= 1,  and  cris  the  standard  parameter 
of  the  primary  Rayleigh  distribution  of  large  mesospheric  electric  fields. 

Using  ft  =  99  samples  of  the  effective  electron  collision  frequency  collected  in  the  height  range  61-67  km 
in  Canada,  the  histogram  showing  the  distribution  of  the  disturbances  in  the  effective  electron  collision 
frequency  is  constructed  (Figure  3  in  [9]).  This  histogram  gives  an  estimate  of  the  77-distribution  first 
ordinary  moment  Mi [77]  =  2.42  ±  0.23  within  the  0.98  confidence  interval  for  the  60  -  67-km  altitude 
range.  These  data  have  been  used  for  inferring  the  electric  fields  whose  distribution  is  presented  in  [9] 
(Figure  6).  The  fitted  theoretical  distributions  j{rj).  Equation  (23),  have  the  following  parameters:  S\  /  0  = 
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0.54  and  ve0  =  3.32  x  107  s  1  for  61  km,  S\  /  o2  =  0.54  and  ve0  =  2.21  x  107  s  1  for  64  km,  Si  /  o2  =  0.18  and 
ve0  =  1.47  x  107  s”1  for  67  km. 

Using  n  —  129  samples  of  the  effective  electron  collision  frequency  collected  at  the  altitudes  of  60-km,  63- 
km,  and  66-km  altitude  over  Kharkiv  (Ukraine),  the  histogram  showing  the  distribution  of  the 
disturbances  in  the  effective  electron  collision  frequency  is  constructed  in  a  similar  fashion  (Figure  4  in 
[9]).  This  histogram  gives  an  estimate  of  the  77-distribution  first  ordinary  moment  Mi[r/\  =  2.02  ±  0.14 
within  the  0.98  confidence  interval  for  the  60  -  66-km  altitude  range.  The  fitted  theoretical  distributions 
ftrj).  Equation  (23),  have  the  following  parameters:  S\  /  a2  =  0.94  and  ve0  =  3.75  x  107  s”1  for  60  km,  Si  /  0 
=  0.84  and  ve0  =  2.55  x  107  s”1  for  63  km,  S)  /  o2  =  0.45  and  ve0  =  1.68  x  107  s”1  for  66  km.  These  data  have 
been  used  for  inferring  the  electric  fields  whose  distribution  is  presented  in  [9]  (Figure  5). 

The  same  datasets  are  used  to  construct  histogram  showing  the  distribution  of  the  disturbances  in  the 
electron  temperature,  we.  The  ^-distribution  first  ordinary  moment  Mx[6\  for  Canada  is  equal  to  2.91  ± 
0.33,  and  2.35  ±0.19  for  Ukraine  [9]. 

The  analysis  of  the  data  on  the  relative  disturbances  in  the  effective  electron  collision  frequencies  and  in 
the  electron  temperatures  shows  that  the  large  mesospheric  electric  fields  from  the  Rayleigh  generator 
maintain  the  electrons  in  the  lower  part  of  the  ionospheric  D  region  at  elevated  temperatures,  a  factor  of  2 
higher  than  Te0  and  the  neutral  temperatures  Tn.  Within  the  0.98  confidence  interval,  the  disturbed  Te  and 
values  at  higher  geomagnetic  latitudes  are,  on  average,  higher  than  at  mid  geomagnetic  latitudes.  The 
information  on  temperatures  is  of  major  importance  in  determining  chemical  reactions  rates. 


3.0  MODELING  DISTURBANCES  IN  CONDUCTIVITY  PROFILES 

The  accuracy  of  maritime  navigation  systems  using  very  low  frequency  signals,  such  as  Omega,  depends 
on  knowing  accurately  the  altitude  of  the  bottom  of  the  ionosphere.  It  is  well  known  that  rapid  vertical 
changes  in  this  boundary  during  solar  flares  and  geomagnetic  storms  can  introduce  errors  of  several 
kilometers  in  location  determinations. 

Since  the  phenomenon  of  large  mesospheric  electric  fields  is  supposed  to  be  local  and  random,  this  leads 
us  to  expect  that  the  surface  of  reflection  is  a  rough  statistical  surface. 

The  large  mesospheric  electric  fields  create  disturbances  in  the  mesospheric  conductivity  that  is  governed 
by  equations  (4)  -  (13).  This  section  briefly  summarizes  numerical  simulations  of  the  conductivity 
disturbances.  The  starting  point  for  modeling  is  the  profile  of  an  electric  field.  The  distribution  of  the  large 
mesospheric  electric  fields,  presented  in  Section  2.4.4,  shows  two  parts,  the  main  body  of  0  <  E  <  2.5  V/m, 
and  the  tail  of  E  >  2.5  V/m.  The  main  body  corresponds  to  a  one-parameter  Rayleigh  probability  density 
function.  The  dataset  for  the  tail  is  too  small  to  draw  statistical  inferences,  and  therefore,  these  electric 
fields  have  to  be  modeled  by  rocket  measurements  presented  in  Figure  1.  Since  the  statistical  analysis 
presented  in  Section  2.4.6  shows  that  the  distribution  function  of  large  mesospheric  electric  field  values 
does  not  exhibit  a  noticeable  diurnal  dependence,  the  same  electric  field  profile  may  be  used  to  model 
daytime  and  nighttime  conditions.  Accordingly,  the  mesospheric  electric  fields  have  been  modeled  by  two 
profiles,  along  with  the  corresponding  electron  density  profiles. 

3.1  Electric  Fields  in  the  0  <  E  <  2.5  V/m  Range 

The  electric  field  profile  representing  the  main  body  of  the  electric  field  distribution  0  <  E  <  2.5  V/m  at  the 
lower  ionospheric  boundary  at  midlatitudes  [10]  is  shown  in  Figure  4.  Here,  the  electric  field  intensities  do 
not  exceed  1  V/m. 
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Figure  5  depicts  the  profile  of  the  low-frequency  conductivity,  cre(z),  disturbed  by  large  mesospheric 
electric  fields  presented  in  Figure  4  under  daytime  conditions  as  compared  to  the  undisturbed  profile, 
creo(z),  shown  as  dashed  line. 

Under  nighttime  conditions,  the  same  electric  field  profile  in  Figure  4  would  produce  the  disturbances  in 
the  low-frequency  conductivity  profile,  cre(z ),  shown  in  Figure  6. 

The  important  feature  to  note  is  that  although  the  electric  field  of  magnitude  below  1  V/m  may  reduce  the 
low-frequency  conductivity  at  the  lower  ionospheric  boundary  by  50%  under  both  daytime  and  nighttime 
conditions,  it  is  unable  to  produce  a  local  minimum  in  the  conductivity  profile  cre(z). 


Figure  4.  The  electric  field  profile  representing  the  main  body  of  the  electric  field 
distribution  0  <  E  <  2.5  V/m  at  the  lower  ionospheric  boundary  at 
midlatitudes,  as  inferred  from  the  Ukrainian  MF  radar  data. 


3.2  Electric  Fields  in  the  E  >  2.5  V/m  Tail 

Figure  7  shows  the  profile  of  the  low-frequency  conductivity,  cre(z),  disturbed  by  the  large  mesospheric 
electric  fields  presented  in  Figure  1  under  daytime  conditions  as  compared  to  the  undisturbed  profile 
creo(z). 


The  same  electric  field  profile  under  nighttime  conditions  would  produce  the  disturbances  in  the  low- 
frequency  conductivity  shown  in  Figure  8. 

The  low-frequency  conductivities  and  their  disturbances  by  the  electric  fields  of  several  volts/meter 
magnitude  for  daytime  and  nighttime  conditions  differ  by  orders  of  magnitude;  however,  their  respective 
magnitudes  not  only  reduce  by  a  factor  of  a  few  times,  but  also  exhibit  a  local  height  minimum. 

3.3  Conclusions 

The  simulations  illustrated  above  demonstrate  that  temporal  variability  in  the  large  mesospheric  electric 
fields  may  produce  spatial  variability  in  the  conductivity  contours.  If  the  large  mesospheric  electric  fields 
exhibit  peak  intensity  above  60  km,  then  a  peak  of  1  V/m  gives  a  local  rise  in  the  conductivity  contour 
heights  by  nearly  6  km,  while  a  4  V/m  peak  intensity  results  in  a  local  10-km  upward  shift.  It  is  important 
to  notice  that  the  disappearance  of  the  large  mesospheric  electric  fields  causes  the  lowering  of  the 
ionospheric  conduction  contours  by  the  same  amount. 
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Conductivity  ( Siemens  per  Meter) 


Figure  5.  The  profile  of  the  low-frequency  conductivity  disturbed  (solid  line)  by  the  large 
mesospheric  electric  fields  with  the  1  V/m  peak  shown  in  Figure  4  under  daytime 
conditions  and  the  undisturbed  (dashed  line)  conductivity  profile. 


Conductivity  {Siemens  per  Meter) 


Figure  6.  The  profile  of  the  low-frequency  conductivity  disturbed  (solid  line)  by  the  large 
mesospheric  electric  fields  with  the  1  V/m  peak  shown  in  Figure  4  under  nighttime 
conditions  and  the  undisturbed  (dashed  line)  conductivity  profile. 
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Conductivity  {Siemens  per  Meter) 


Figure  7.  The  profile  of  the  low-frequency  conductivity  disturbed  (solid  line)  by  the  large 
mesospheric  electric  fields  with  the  4  V/m  peak  shown  in  Figure  1  and  the  undisturbed 
(dashed  line)  conductivity  profile  under  daytime  conditions. 


SS 


Conductivity  {Siemens  per  Meter) 


Figure  8.  The  profile  of  the  low-frequency  conductivity  disturbed  (solid  line)  by  the  large 
mesospheric  electric  fields  with  the  4  V/m  peak  shown  in  Figure  1  and  the  undisturbed 
(dashed  line)  conductivity  profile  under  nighttime  conditions. 
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4.0  VLF  PHASE  PERTURBATIONS  PRODUCED  BY  LARGE  MESOSPHERIC 
ELECTRIC  FIELDS 

The  calculations  of  VLF  phase  variations  require  3D  models  of  the  ionospheric  conduction  contours, 
which  remain  unknown.  Therefore,  this  section  is  concerned  with  the  most  logical  ways  to  proceed. 


4.1  Smooth  Reflective  Surface 

Since  the  data  collected  thus  far  provide  only  a  one-dimensional  view  of  the  vertical  variability  in  the  level 
of  reflection,  the  VLF  phase  perturbations  could  be  assessed  correctly  only  if  the  reflective  surface  is 
smooth.  If  the  mean  height  of  reflection  h  changes  by  an  amount  Ah,  the  corresponding  phase  change  A</> 
over  long  distances  can  be  given  (in  radians)  by  [36] 


A  (j>~- 


2nd 


X 


h  X 
-  + 


2  A 


2  a  16  r 


Ah 


(25) 


where  d  is  the  great  circle  distance  of  the  receiver  from  the  transmitter,  and  a  is  the  radius  of  the  Earth. 
The  wavelength,  X,  of  special  interest  is  approximately  equal  to  30  km,  while  vertical  excursions  in  the 
contours  of  constant  conductivity  due  to  the  large  mesospheric  electric  fields  are  estimated  to  be  less  than 
10  km.  This  way  of  thinking  was  employed  in  interpreting  the  VLF  phase  variations  (at  16  kHz)  attributed 
to  variations  in  the  large  mesospheric  electric  fields  during  the  escape  of  radioactive  materials  [37]. 


4.2  Statistically  Rough  Reflective  Surface 

If  the  contours  of  constant  conductivity  that  determine  the  phase  fluctuations  of  the  reflected  VLF  signals 
are  statistically  rough  surfaces,  then  the  determination  of  the  fluctuations  of  the  signal  parameters  require 
numerical  information  on  the  character  of  surface  irregularity.  The  theory  of  wave  scattering  from 
statistically  rough  surfaces  requires  information  at  least  about  the  correlation  functions  of  the  deviations  of 
reflective  surface  from  the  smooth  surface  [38].  It  is  therefore  evident  that  in  order  to  adequately 
determine  the  effects  of  large  mesospheric  electric  fields  on  radiowave  propagation,  the  variations  in  the 
surfaces  of  constant  VLF  reflectivity  must  be  studied  and  specified  much  more  accurately. 

However,  the  MF  radars  used  so  far  do  not  permit  the  determination  of  the  three-dimensional  distribution 
of  electric  fields  and  do  not  provide  continuous  real-time  observations  of  large  mesospheric  electric  fields. 
The  first  restriction  is  due  to  the  MF  radar  antennas  forming  wide  fixed  vertical  beams  and  thus  capable  of 
measuring  only  a  one-dimensional  distribution  of  ionospheric  parameters.  The  second  restriction  is  due  to 
the  low  signal-to-noise  ratios,  and  consequently,  a  new  MF  radar  facility  should  provide  higher  signal-to- 
noise  ratios. 

Thus,  a  challenge  for  VLF  phase  calculations  will  be  to  characterize  three-dimensional  variations  in  the 
electric  fields  intrinsic  to  the  mesosphere.  Consequently,  this  task  presents  a  challenge  for  the  construction 
of  an  MF  radar  facility  capable  of  studying  the  three-dimensional  nature  of  large  mesospheric  electric 
fields. 

The  engineering  solution  to  the  problem  of  determining  a  3-dimentianal  distribution  of  scattered  signals  in 
the  MF  frequency  band  has  yet  to  be  found.  At  least  two  ways  of  obtaining  this  solution  may  be  indicated. 
First,  a  3-dimentianal  distribution  could  be  provided  by  a  new  radar  capable  of  performing  elevation  and 
azimuth  scans.  A  simpler,  but  not  a  necessarily  better,  solution  may  be  a  few  relocatable  MF  radars  that 
may  be  relocated  in  response  to  new  knowledge  and  understanding  of  the  determining  factors  underlying 
mesospheric  electrodynamics. 
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5.0  CLUSTERED  INSTRUMENT  STUDIES 

Since  the  construction  of  the  new  MF  radar  facility  requires  funds,  which  are  not  available  at  present,  joint 
co-located  MF  radar  and  VLF  propagation  measurement  campaigns  may  provide  some  information  useful 
for  practical  applications.  It  is  evident  that  one-hop  VLF  circuits  with  an  MF  radar  below  the  region  where 
the  ray  is  reflected  may  yield  new  advances  in  VLF  phase  perturbation  research.  The  Kharkiv  V.  N. 
Karazin  National  University  developed,  produced,  and  used  relocatable  MF  radars  in  numerous  clustered 
instrument  measurement  campaigns  at  high  and  middle  latitudes  in  the  U.S.S.R.  [17]. 

However,  the  electrodynamics  of  the  mesosphere  is  one  of  the  most  poorly  understood  topics,  and  its 
studies  involve  mesospheric  winds  and  complicated  chemistry  of  the  mesosphere  and  of  the  ionospheric  D 
region.  These  studies  require  coordinated  efforts  of  the  entire  scientific  community  who  could  utilize  all 
radio  and  optical  techniques  and  computer  simulations  to  cover  a  wide  variety  of  ionospheric  conditions. 


6.0  CONCLUSION 

The  analysis  of  the  complex  nature  of  the  problems  encountered  in  assessing  VLF  phase  perturbations 
produced  by  the  large  mesospheric  electric  fields  has  provided  the  following  inferences. 

The  large  mesospheric  electric  fields  act  to  produce  local  height  variations  in  the  ionospheric  conduction 
contours  of  the  order  of  a  few  kilometers. 

A  major  achievement  in  this  area  has  been  the  development  of  the  MF  radar  technique  for  sensing  the 
large  mesospheric  electric  fields  remotely. 

The  MF  radars  used  so  far  have  two  restrictions.  The  minor  restriction  is  low  levels  of  the  signal-to-noise 
ratio  in  the  60-70  km  altitude  range,  particularly  under  nighttime  conditions,  whereas  the  technique 
requires  the  signal-to-noise  ratio  to  be  equal  or  greater  than  about  five.  An  increase  in  the  signal-to-noise 
ratio  is  the  way  to  continuously  monitor  large  mesospheric  electric  fields. 

The  fundamental  restriction  is  the  fact  that  the  existing  MF  radars  provide  only  one-dimensional 
distribution  of  the  electric  fields,  while  the  modeling  of  VLF  phase  perturbations  requires  information  on 
the  3-dimensional  variability  of  the  ionospheric  conduction  contours.  The  absence  of  the  3-dimensional 
distribution  of  the  large  mesospheric  electric  fields  is  the  major  obstacle  to  progress  in  studying  VLF 
phase  perturbations.  Therefore,  the  most  important  challenge  for  the  future  will  be  to  develop  an  MF  radar 
facility  capable  of  determining  a  3-dimentianal  distribution  of  scattered  signals  in  the  MF  frequency  band. 
The  data  this  facility  will  acquire  in  some  area  of  the  Earth  would  be  sufficient  to  model  VLF  phase 
perturbations  produced  by  the  large  mesospheric  electric  fields  in  that  area. 

Eventually,  an  overall  understanding  of  the  dynamics  and  mesospheric  and  ionospheric  D-region 
chemistry,  which  establish  conductivity  patterns,  can  be  achieved  by  the  combined  efforts  of  the  entire 
scientific  community. 

Meanwhile,  one-hop  VLF  circuits  with  an  MF  radar  below  the  region  where  the  ray  is  reflected  may  yield 
new  advances  in  VLF  phase  perturbation  research  and  provide  some  information  useful  for  practical 
applications. 
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Historical  Review 

First  Measurements 
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Rocket-borne  field  mill. 

Yu.  A.  Bragin,  A.A.  Tyutin, 
A.A.  Kocheev,  and  A.A.  Tyutin 
Direct  measurement  of  the 
atmospheric  vertical  electric 
field  intensity  up  to  80  km. 
Cosmic  Res.,  12.  279,  1974. 


Historical  Review 

Rocket  Electric  Field  Data 


N.C.  Maynard,  L.C.  Hale, 

J.D.  Mitchell,  F.J.  Schmidlin, 
R.A.  Goldberg,  J.R.  Barcus, 
F.  Soraas,  and  C.L.  Croskey, 
Electrical  structure  in  the 
high  latitude  middle 
atmosphere, 

J.  Atmos.  Terr.  Phys.,  46, 
NqI,  1984. 


Historical  Review 

Worldwide  Rocket  Database 


1 974  -  2006 

Less  than  100  rockets  launches 

Conclusion 

Not  enough  launch  data  to  address  VLF 
phase  perturbations 

To  achieve  a  breakthrough,  a  remote  sensing 
technique  employing  radio  waves  is  required 


Remote  Sensing  Technique 

Instrumentation 

MF  Radar: 

2.2  MHz  operating  frequency 
20-ps  pulse 

300  x  300  m2  transmitter  antenna 
Receiving  array  of  two  circularly  polarized 
antennas 


Remote  Sensing  Technique 

Data  Acquired 

Magnetoionic  Components: 

•  Amplitudes  of  circularly 

polarized  waves 

•  Amplitudes  of  circularly  polarized 

waves 

•  Backscattered  from  60-  to  70-km  altitude 
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Remote  Sensing  Technique 

Data  Input 


The  input  data  for  our  MF  radar  technique 
consists  of  the  ratios  of  the  sample  ordinary 
and  extraordinary  amplitudes  squared 

Backscattered  from  60-  to  70-km  altitude 


Remote  Sensing  Technique 

MF  Radars  Used  at  Present 

Measured  vs.  Required  Features 

Signal-to-noise  ratio  regularly  falls  below  the 
acceptable  level 

•  Temporal  continuity 

One-dimensional  distribution 


3D  distribution 


Remote  Sensing  Technique 

Governing  Relations 

Weakly  Ionized  Plasma: 

Electrons,  positive  and  negative  ions,  neutrals 

•  Electron  and  ion  continuity  equations 

•  Energy  equation 

•  Ohm's  law  for  electric  fields 
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Remote  Sensing  Technique 

Basic  Relation  Derived 


•  We  have  related  these  ratios  to  the  electric 
field  intensity  via  the  electron  collision 
frequency 

•  We  have  obtained  the  following  relation  for 
the  electric  field  intensity: 


Remote  Sensing  Technique 

Electron  Collision  Frequency 

First,  we  apply  the  “Differential  Absorption” 
technique  of  [Belrose,  J.  S.,  and  M.  J. 
Burke,  A  study  of  the  lower  ionosphere 
using  partial  reflections.  1.  Experimental 
technique  and  method  of  analysis,  J. 
Geophys.  Res., 

Vol.  69,  No.  13,  pp.  2799-2818,  1964]  at 

altitudes  of  60  to  70  km  and  determine  the 
electron  collision  frequencies 


Remote  Sensing  Technique 

MF  Radar  Electric  Field  Database 

Canada,  1979-1982 

measurements,  duration 

Ukraine,  1978  - 1994 

measurements,  duration 


Total  MF  radar  database  of 
355  measurements 


Remote  Sensing  Technique 

MF  Radar  Measurements 

Temporal  Variability 

The  database  consists  of  separate  intervals 
of  measurements  collected  over  the  years, 
therefore  the  temporal  variability  they 
provide  is  not  continuous 

Nevertheless,  brief  snapshots  can  be 
provided 


Remote  Sensing  Technique 

MF  Radar  Measurements 


Temporal 

•  The  variability  of  large 
mesospheric  electric 
fields  displays  a  wide 
range  of  time  scales, 
from  an  order  of  one 
minute  (as  shown  in 
this  figure)  to  ... 


Variability 


15 


Remote  Sensing  Technique 

MF  Radar  Measurements 

Temporal  Variability 


...an  order  of  an 
hour  (as  shown  in 
this  figure) 


Remote  Sensing  Technique 

MF  Radar  Measurements 

Distribution  of  Electric  Field  Intensities 

Two  parts: 

•  the  main  body  of 
0  <  E  <  2.5  V/m 

(n  =  99) 

•  the ‘tail’ E>  2.5  V/m 

(n  =  21) 


E  {V/m) 
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Remote  Sensing  Technique 

MF  Radar  Measurements 

Distribution  of  Electric  Field  Intensities 


Two  parts: 

•  the  main  body  of 
0  <  E  <  2.5  V/m 

(n  =  129) 

•  the ‘tail’ E>  2.5  V/m 

(n  =  10) 


E  (V/m) 
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Remote  Sensing  Technique 

MF  Radar  Measurements 

Distribution  of  Electric  Field  Intensities 


•  The  main  body 
corresponds  to  a 
one-parameter 
Rayleigh  probability 
density  functions 
as  given  by 
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Disturbances  in  Conductivity 

Electric  Field  Model 


To  model 
disturbances  in 
conductivity  profiles 
produced  by 

fields,  we  have 
chosen  the  profile 
depicted  in  this 
figure 


Disturbances 

(Daytime 

•  Shown  here  is 
the  profile  of  the 
low-frequency 
conductivity 
disturbed  (solid  line) 
by  the 

electric  fields 
and  the  undisturbed 
(dashed  line) 
conductivity  profile 


in  Conductivity 

Conditions) 
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Disturbances  in  Conductivity 

(Nighttime  Conditions) 


Shown  here  is 
the  profile  of  the 
low-frequency 
conductivity 
disturbed  (solid  line) 
by  the  less  than  1 
electric  fields 
and  the  undisturbed 
(dashed  line) 
conductivity  profile 


iO"10  i  cr1'  iO“8  10“7  lo-6 

Conductivity  ( Siemens  per  Meter ) 


Disturbances  in  Conductivity 

Conclusion 


•  The  electric  field  variability  in  the  main  body 
of  the  distribution  leads 

to  variations  in  the  ionospheric  conduction 
contours  of  up  to  a  few  kilometers 

in  altitude 
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Disturbances  in  Conductivity 

Electric  Field  Model 


To  model 
disturbances  in 
conductivity  profiles 
produced  by 

fields,  we  have 
chosen  the  mean  of 
the  profiles  depicted 
in  this  figure 
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Disturbances  in  Conductivity 

(Daytime  Conditions) 


Shown  here  is  the 
profile  of  the  low- 
frequency  conductivity 
disturbed  (solid  line)  by 
the  E  >  2.5  V/m  electric 
fields  shown  in  the 
preceding  figure  and 
the  undisturbed 
(dashed  line) 
conductivity  profile. 
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Disturbances  in  Conductivity 

(Nighttime  Conditions) 

Shown  here  is  the 
profile  of  the 
low-frequency 
conductivity 
disturbed  (solid  line) 
by  the 

electric  fields  and 
the  undisturbed 
(dashed  line) 
conductivity  profile 


Disturbances  in  Conductivity 

Iso-Conduction  Contours 

Conclusions 

The  disappearance  or  appearance  of  electric 
fields  in  the  main  body  of  their  distribution  (with 
a  peak  of  1  V/m)  results  in  a  variation  in  the 
height  of  the  ionospheric  iso-conduction 
contours  of  up  to  6  kms  in  altitude 

The  disappearance  or  appearance  of  electric 
fields  in  the  “tail”  of  their  distribution,  with  a  peak 
of  4  V/m,  results  in  a  variation  in  the  height  of 
the  ionospheric  iso-conduction  contours  of  up 
to  10  kms  in  altitude 


VLF  Phase  Perturbations 

(Smooth  Reflective  Surface) 


If  the  reflective  surface 
were  smooth,  the  VLF 
phase  perturbations 
would  be  assessed 
correctly.  The 
corresponding  phase 
change  over  long 
distances  can  be  given 
(in  radians)  by  the 
well-known  relation 


2nd  [  h  A2  )Ah 
/l  i  2 a  16/r  J  h 


VLF  Phase  Perturbations 

(Smooth  Reflective  Surface) 

•  Given  this  way  of  thinking,  we  have  attributed 
the  VLF  phase  variations  (at  16  kHz)  to 
variations  in  the  electric  fields  during  the 
escape  of  radioactive  materials  I.M.  Fuks, 
R.S.  Shubova,  and  S.l.  Martynenko, 

Lower  ionosphere  response  to 
conductivity  variations  of  the  near-earth 
atmosphere, 

J.  Atmos.  Solar-Terr.  Phys.,  Vol.  59,  961, 
1997] 
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VLF  Phase  Perturbations 

(Rough  Reflective  Surface) 

•  The  theory  of  wave  scattering  requires 
information  on  the  correlation  functions  of  the 
deviations  of  the  reflective  surface  from  a 
smooth  surface 

•  Consequently,  the  challenge  for  VLF  phase 
investigations  will  be  to  characterize  three- 
dimensional  variations  in  the  electric  fields 
intrinsic  to  the  mesosphere 
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VLF  Phase  Perturbations 

(Rough  Reflective  Surface) 

Engineering  Solutions 

•  A  new  radar  capable  of  performing  elevation 
and  azimuth  scans  (too  complicated!) 

•  A  few  relocatable  MF  radars  (ones 
developed,  produced  and  used  by 
Kharkiv  V.  Karazin  National  University) 
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Clustered  Instrument  Studies 


•  One-hop  VLF  circuits  with  an  MF  radar 
below  the  region  where  the  ray  is  reflected 
may  yield  new  advances  in  VLF  phase 
perturbation  research 
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Conclusions 


The  mesospheric  electric  fields  act  to 
produce  local  height  variations  in  the 
ionospheric  conduction  contours  of  the 
order  of  a  few  kilometers 

A  major  breakthrough  in  this  area  has  been 
the  development  of  the  MF  radar  technique 
for  sensing  electric  fields  remotely 


Conclusions 


The  MF  radars  used  so  far  have  low  levels  of 
the  signal-to-noise  ratios,  which  need  to  be 
increased 

The  existing  MF  radars  provide  only  one¬ 
dimensional  distribution  of  the  electric  fields. 
Thus,  the  most  important  challenge  for  the 
future  will  be  to  develop  an  MF  radar  facility 
capable  of  determining  a  3-dimensional 
distribution  of  scattered  signals  in  the  MF 
frequency  band 


Conclusions 


•  One-hop  VLF  circuits  with  an  MF  radar  below 
the  region  where  the  ray  is  reflected  may 
yield  new  advances  in  VLF  phase 
perturbation  research  and  provide  useful 
information  for  practical  applications 
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